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ABSTRACT 

We perform self-consistent Hartree calculations for a model of A^ x Cqq (A= K, 
Rb) to study the effects of the presence of vacancies. We find that the strong vacancy 
potential is very efficiently screened and that the density of states is only weakly 
influenced by the presence of vacancies. 

I. INTRODUCTION 

The A3C60 (A= K, Rb) compounds with the nominal stoichiometry three are 
believed to actually have the composition A3_ x .Ceo, with x ~ . T . The alkali 
atoms are believed to be almost completely ionized, and the electrostatic potential 
from the alkali ions should therefore be large. The presence of a vacancy may then 
be a large perturbation on the electronic structure. For instance, if a vacancy were 
completely unscreened, its potential at the centre of a neighboring C 60 molecule would 
be about 2.3 eV and it would split the t\ u states by about 0.2 eV, which is comparable 
to the t\ u band width (~ 1/2 eV) for a system without vacancies. The vacancies occur 
on the tetrahedral positions |TJ and each vacancy is surrounded by four Ceo molecules. 
If x = 0.07, more than a quarter of the Ceo molecules have an alkali nearest neighbor 
vacancy. It is an interesting question how this influences the electronic structure||. 
Similar questions arise for the systems Na2Cs x C 60 and Ka.^Ba^CeotQ, which also 
have different charges at the alkali positions. 

We perform Hartree calculations for a model of As^Ceo- The model includes 
the 60 "radial" 2p orbitals on each Ceo molecule, which are allowed to screen the 
potential from the alkali ions self-consistently. We find that the screening is very 
efficient, in particular in terms of charge transfer between the molecules and to a less 
extent in terms of polarizing the molecules. As a result the density of states is almost 
unchanged by the presence of the vacancies. 

In Sec. II we present the model, in Sec. Ill the method for performing the 
calculations. The results are given in Sec. IV and some implications are discussed in 
Sec. V. 



II. MODEL 

We introduce a model of As-^Ceo for which we calculate the electronic structure 
and the screening properties. For each carbon atom in each molecule we include the 
2p orbital, pointing radially out from the molecule. The states close to the Fermi 
energy are mainly formed from these radial 2p orbitals, and they are therefore par- 
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ticularly important for the screening. The cr-like 2p orbitals pointing tangentially to 
the molecular surface and the 2s orbitals contribute mainly to orbitals far away from 
the Fermi energy and are therefore less important. If these orbitals had nevertheless 
been included, the result would have been an even more efficient screening than found 
below. Thus we consider the one-particle Hamiltonian 

ivu ivj/MT 

where e gives the energy of the 2p states and t(iu,jfi) gives the hopping integrals 
between these states. The atoms in a Cqo molecule are labelled by Roman letters and 
the C 6 o molecules are labelled by Greek letters. This includes both hopping within a 
molecule and hopping between states on nearest neighbor molecules. We have used 
the parametrization in Ref. but multiplied all hopping integrals by a factor 1.2 to 
obtain a t\ u band width in agreement with recent band structure calculations. The 
Ceo molecules have an orientational disorder, and take essentially randomly one out 
of two preferred directions f|. This is built into the hopping matrix elements ||. 
The Coulomb integrals between the orbitals are defined as 

rfVrfV0(r-R^)0(r-R jM )|— — 7 |0(r'-R m7 )0(r'-R n<5 ), (2) 

where H iu is the position of the zth atom on the uth C 60 molecule. We neglect the 
overlap of the functions centred on different atoms, and therefore v is nonzero only if 
(iu) = (jfi) and [m^) = (n8). We make the assumption 

for R;„ — Rm/J > 



v(iv, jfi;mj,n8) = S itj 5 u>fl S mtn 8 7tS * I l R - R ^l " m7 (3) 

{ v for \H iu - R m7 | = 

For the on-site interaction, we have obtained vq = 12 eV from atomic calculations. |7j 
For the Coulomb integrals between an alkali core and the carbon 2p orbitals a similar 
expression is used. 

The Ceo molecules are put on a fee lattice with the lattice parameter 14.24 A, 
and the alkali atoms are placed on the octahedral and tetrahedral positions. It was 
assumed that all Cgo molecules and alkali ions are at their ideal positions, and that 
there are no distortions due to the vacancies. We introduce a large supercell with 500 
molecules, which is periodically repeated. This supercell is large enough to give results 
which are essentially converged with respect to the cell size. Within the supercell 
each Cgo molecules is allowed to take one of the two preferred orientations randomly. 
Furthermore vacancies are introduced randomly on the tetrahedral positions in such 
a way that on the average there alkali atoms per unit cell. To keep the 

system neutral, we fill up the t\ u band with exactly as many electrons as there are 
alkali ions in the system. 
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III. METHOD 



We perform self-consistent Hartree calculations. For strongly correlated systems, 
this approach might be expected to give poor description of the screening. By com- 
paring Quantum Monte Carlo and RPA calculations, it has, however, been found 
that RPA remains surprisingly accurate up values of the Coulomb interaction where 
the Mott-Hubbard transition takes placed], ||. Since even stoichiometric A3C60 is 
believed to be on the metallic side of a Mott-Hubbard transition 0, we expect the 
Hartree approximation to be qualitatively correct in the present situation. 

To study the model described above, we use the molecular solid character of 
the system. The potential from the vacancies is screened in two ways. First, there 
is a charge transfer between the molecules to approximately neutralize the alkali 
vacancies, i.e., molecules which are nearest neighbors to a vacancy have a smaller 
electronic charge. This screening is very efficient and mainly takes place in the partly 
filled t\ u band. Below we only consider the charge transfer in this band. Second, 
the molecules polarize so that within a given molecule the charge moves away from 
a vacancy. Below we assume that the polarization can be calculated by studying an 
isolated molecule in the potential of the surrounding molecules and alkali ions. 

First we calculate the potential for each atom in each molecule due to the the 
alkali ions and the charges on all other atoms in the Cgo molecules. In this process 
we take into account that the "core" of a carbon atom has the charge +1, since the 
a electrons are counted in the core. This calculation is performed using a Madelung 
summation technique. For each molecule we then solve the corresponding 60 x 60 
Hamiltonian, to obtain the 7r-orbitals. The lowest 30 orbitals are occupied, and the 
corresponding charges are calculated for the 60 atoms in each molecule. This describes 
the polarization of the molecules. 

To describe the charge transfer between the molecules, we consider the three t\ u 
orbitals on each atom. From the molecular calculations described above, we obtain the 
energies of the t\ u orbitals, which in general are not degenerate due to the random 
arrangement of the alkali vacancies. These energies give the diagonal elements in 
a 3iV x 3N matrix, where iV is the number of molecules in the unit cell. From 
the calculated coefficients of the t\ u orbitals and the hopping integrals between 2p 
orbitals on different molecules, we obtain the hopping matrix elements between t\ u 
orbitals on different molecules. These give the nondiagonal matrix elements in the 
3iV x 3N Hamiltonian matrix. We use periodic boundary conditions when setting up 
the hopping part of the matrix, which corresponds to considering the k = solution. 
Due to the large size of the supercell and the corresponding small size of the Brillouin 
zone, this should be a good approximation. Finally, we diagonalize the 3N x 3iV 
matrix and obtain the t\ u contribution to the charge density. 
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Figure 1: Density of states (per eV-mol-spin) for A 3 _ x C6o for x = (full curve) and 
x = 0.07 (dashed curve). 

IV. RESULTS 

The results for the density of states N(e) are shown in Fig. |l|. Some small wiggles 
in N(e) are due to the finite cell size (500 molecules). This does not, however, influence 
the conclusions. The main result of Fig. [I] is that N(e) is practically identical for 
x = and x = 0.07. This small difference is due to the efficient screening of the alkali 
vacancies. 

To discuss the screening we introduce the average position 

1 3 

= - Z) e *i«( im )> ( 4 ) 

6 rn=l 

and standard deviation 

1 3 

otu. W 2 = o J2 ( £ hu (im) - e tlu (i)) 2 , (5) 

6 m=l 

of the t\ u level, where Et lu (im) is the energy of the mth orbital on the ith molecule. 
We also introduce the average of this standard deviation 

j N 

a tiu = T7 5Z (7 ( i )- ( 6 ) 

ly i=l 

cr(i) tlu and cr tlu measure how the degeneracy of the t lu orbitals is lifted due to the 
random positions of the vacancies. It is also interesting to see how the average 
ti u position varies between the sites. For this purpose we introduce the standard 
deviation 

1 N 

K, = TrY,(^(i)-etJ 2 , (7) 

1=1 

where Et lu is the average of the t\ u energies over the super cell. In a similar way 
we define A tetra and A oct as the standard deviation in the potentials at the occupied 
tetrahedra and octahedral positions, respectively. 
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Table 1: On-site splitting a tlu (Eq. (^)) and off-site splitting A tlu (Eq. (^)) of the t\ u 
level as well as the splitting in the potential on the tetrahedral (A tefra ) and octahedral 
(A oct ) occupied positions for A 3 _ a .C 60 . The first three lines refer to the system without 
vacancies (x = 0) and the following three lines to a system with 7% (x = 0.07). The 
band width W is also shown. All energies are in eV. 
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A tlu 
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Nonself-consistent 


0.00 


0.00 


0.00 


0.00 


0.00 


0.53 


t\ u self-consistent 


0.0 


0.00 


0.01 


0.01 


0.01 


0.54 


Self-consistent 


0.00 


0.00 


0.01 


0.01 


0.01 


0.55 


Nonself-consistent 


0.07 


0.31 


6.83 


6.76 


6.72 


32.31 


t\ u self-consistent 


0.07 


0.01 


0.04 


0.08 


0.20 


0.71 


Self-consistent 


0.07 


0.01 


0.02 


0.01 


0.05 


0.55 



As discussed above, there are two types of screening. Due to rearrangements 
of the partly occupied t\ u band, charge is transferred between the molecules. This 
influences the averaged positions e tlu (i) of the t\ u level on the different molecules, 
and it therefore reduces A tlu . This also reduces the splitting of the t\ u levels inside 
a molecule. This splitting is further reduced by rearrangements of the charge inside 
the molecules. 

We therefore perform three types of calculations. First, we perform a nonself- 
consistent calculation, where only the potentials from the alkali ions are considered, 
and any screening from the C 60 molecules is neglected. Second, we perform a cal- 
culation where charge transfer between the molecules is allowed, but polarization 
of the Ceo molecules is neglected (t\ u self-consistent). Finally, we perform the full 
self-consistent calculation as discussed above. The results are summarized in Table 

i 

In the nonself-consistent calculation, we find a very large variation (A tlu ) in the 
average position of the ti u level between different molecules for x = 0.07. There is 
also a large splitting (cr tlu ) of the t\ u level on molecules close to a vacancy. As a result 
the band width W is huge. When charge transfer between the molecules is allowed 
self-consistency) the variation in the t\ u level position between the molecules 
becomes very small. In the fully self-consistent calculation, the polarization of the 
molecules reduces the splitting of the ti u level somewhat more. Thus the largest value 
of <Jtiu{i) is reduced by about a factor of four and the average value by about a factor 
of two. Due to the rounding in Table ^ this is not seen explicitly. Similar results are 
found for the potential on the alkali ion positions, and the potential on the different 
occupied tetrahedral sites only varies by a few hundredths of an eV. 
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V. CONCLUDING REMARKS 



We find that the strong potential from an alkali vacancy is efficiently screened 
by charge transfer between the molecules and to some extent also by polarization 
of the molecules. As a result, the density of states N(e) is hardly influenced by 
the vacancies. These results were obtained under the assumption that there are no 
distortions of the Ceo molecules and alkali ions around the vacancies. 

One might have expected that the strong vacancy potential would lift the t\ u 
degeneracy on a given molecule, thereby increasing the band width with vacancy 
concentration. This would have reduce N(0) and the electron-phonon interaction A, 
and it might have explained the rapid drop|| in T c as the stoichiometry deviates 
from three. Due to the efficient screening found above, such an explanation, however, 
seems unlikely. 

In NMR of Rb3_a;C6o (x ~ 0.03) it is found that the quadrupolar distortion of the 
Rb lines is rather small, in spite of the electric field gradients one would expect from 
the vacancies 0]. It was argued that this may be due due to an efficient screening of 



the vacancies [10]. Such considerations agree very well with the present results. 
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